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FOREWORD 


This  interim  technical  report  was  prepared  by  the  Yale  University  Observa¬ 
tory,  New  Haven,  Connecticut  on  Contract  AF  33(657)-7271  for  the  Aeronautical 
Research  Laboratories,  Office  of  Aerospace  Research,  United  States  Air  Force. 
The  research  reported  herein  was  accomplished  on  Task  707  3-01,  "High  Energy 
Generation  and  Control"  of  Project  7073,  "Research  on  Plasma  Dynamics  under 
the  Technical  Cognizance  of  Dr.  W.  G.  Braun  of  the  Plasma  Physics  Research 
Laboratory  of  ARL. 


11 


ABSTRACT 


A  general  formulation  of  the  theory  of  departures  from  LTE  in  free-free 
emissions  and  absorption  is  developed.  The  problem  is  treated  subject  to  the 
restrictions  to  steady  state  and  radiative  interactions.  Particle  distribution 
function  and  radiation  field  are  studies  separately,  special  attention  is  given  to 
the  transfer  problem.  The  possibility  of  negative  absorption  (maser  action)  is 
discussed. 
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1,  INTRODUCTION 


The  importance  of  departures  from  local  thermodynamic  equilibrium  (LTE) 

was  first  recognized  in  connection  with  certain  astrophysical  problems,  in 

particular,  in  connection  with  the  interpretation  of  the  spectrum  of  the  outer  solar 

atmosphere \  Ftellar  plasmas  are,  in  the  observable  layers,  rather  weakly 

ionized,  so  that  the  major  part  of  the  radiation  spectrum  consists  of  spectral 

lines,  i,  e,  ,  emissions  and  absorptions  involving  only  bound  states,  with  the  free 

energy  states  in  LTE.  In  physical  terms  this  statement  implies  that  the  free 

energy  states  of  the  electrons  are  distributed  according  to  a  Maxwellian  distribu~ 

2 

tion,  whereas  emission  and  absorption  are  connected  by  Kirchhoff’s  law  , 

In  many  domaines  of  laboratory  physics,  however,  plasmas  are  completely 
ionized,  and  the  population  of  bound  states  can  be  altogether  neglected.  Then, 
only  the  free  electrons  are  involved  in  radiative  interactions,  and  departures 
from  a  Maxwellian  distribution  of  the  free  energy  states. 

The  purpose  of  this  paper  is  to  outline  a  general  formulation  of  departures 
from  LTE  in  free-free  emissions  and  absorptions.  In  this  report,  we  restrict 
our  considerations  to  steady  states,  where  the  change  in  the  population  of  states 
is  (during  observation  times)  small  compared  with  the  populations  themselves, 
and  to  radiative  transitions  which  involve  emission  or  absorption  of  a  single 
photon. 

Radiative  interactions  affect  the  particle  distribution  and  the  spectrum  of  the 
radiation  field.  There  are  therefore  always  two  sets  of  conditions  to  be  con¬ 
sidered,  derived  from  the  behavior  of  the  radiation  field,  and  from  the  behavior 
of  the  particles.  Neglecting  one  set  of  conditions  does  mathematically  not  lead 
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to  any  contradiction,  but  leads  to  physically  invalid  solutions.  This  fact  has  been 

3 

pointed  out  previously  . 

We  begin  by  deriving  the  steady  state  condition  for  the  particle  distribution 
function  (Section  2),  turning  to  the  radiation  field  in  Section  3,  In  Section  4,  the 
trivial  solution  (LTE)  is  discussed.  Section  5  deals  with  some  obvious  physical 
consequences  of  the  steady* state  condition.  Section  6  contains  a  discussion  of  the 
relations  between  particle  distribution  function  and  local  intensity  of  the  radiation 
field.  In  an  Appendix,  the  question  of  possible  maser  action  involving  continuous 
transitions  is  discussed. 

2.  STEADY-STATE  CONDITION  FOR  THE  ENERGY  DISTRIBUTION  OF  FREE 
PARTICLES 

If  a  steady  state  has  been  reached  in  a  plasma,  the  number  of  particles 
entering  a  certain  energy  state  must  equal  the  number  of  particles  leaving  that 
state.  This  relation  then  holds  for  all  energy  states. 

Consider  an  arbitrary  state  with  energy  E^,  populated  by  N(E^) 
particles  per  unit  volume.  Neglecting  non-radiative  transitions,  the  following 
processes  will  occur:  the  particle  may  absorb  a  quantum 

o 

Si  =  nw  ( 1 ) 

from  the  indident  radiation  field  and  move  up  to  the  energy  level 

E  =  E  +  n, 

1  o 

or  it  may  emit  a  quantum  SI,  descending  to  the  energy  level 

E,  =  E  -Si. 

1  o 


(2) 


(3) 
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The  number  of  these  processes  per  unit  time  depends  on  the  interaction 

mechanisms  (bremsstrahlung,  cyclotron  radiation,  Cerenkov  radiation,  etc.), 

on  the  local  intensity  of  the  radiation  field  (absorption  and  induced  emission  only), 

and  on  the  number  of  particles  N(E)  available. 

The  pertinent  quantities  were  defined  and  discussed  in  detail  in  a  previous 
2 

publication  .  We  recall  that  the  number  of  particles  in  the  energy  range  E,E+dE 
has  the  dimension  [N{E)  dE]  =  cm  so  that  the  LTE  -  distribution  (Maxwell 
distribution)  reads 

N(E)dE=4TTN  ?  e^dE  .  (4) 

o 


Writing  the  radiation  field  too  in  terms  of  energies  ( Ji)  instead  of  angular  frequen¬ 
cies  (w),  we  have  [I_dn]  =  erg/c  m^  sec  sterad,  so  that  in  LTE  (Planckian  distri- 
bution) 


I  =  B^dQ  : 


A  ^ 

4  TT  n  C 


r  «/KT  ,,-l 
[e  -  1]  dfl. 


(5) 


Later  on,  the  abbreviation 


H 

o 


.3  3  2  .  3  2 

4Tr  n  c  [erg  cm 


sec  sterad] 


(6) 


will  be  used. 

The  relations  between  the  cross  sections  of  spontaneous  emission,  absorption 
and  induced  emission  read^ 


H 

o 


a" 


Q 


sp 
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and 


In  the  following, 

WE, a)  s  , 


(9) 


2  .  is  the  energy  of  the  upper  state  from  which  the  spontaneous 
where  E  =  m  v  /2 

emission  takes  place.  The  probability  coefficient  Q  then  has  such  a  dimension 
that  [Q(E,n)df^]=  sec  ^  sterad  \ 

We  are  now  ready  to  write  down  the  number  of  incoming  (n^)  and  outgoing 
(n  )  processes  for  the  arbitrary  energy  state  E^: 


7  H 

=  I  N(Eo  +11)  Q(Eo  +Il,a)  [1  +  ^  I^i]  dO-  + 


+  \  N(Eo 


(10) 


A)  Q(Eo,n)i 


and 


n 


oo 


5  N(E„)  Q(E„+n,n)j3^^ij^dn  + 

+  J  N(E^)  Q(Eo,n)  [1  +  da 


(11) 
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In  steady  state  > 


n 


=  n. 


(12) 


or,  after  some  algebra,  ^ 


^  Q(Eo+n,n.)[Eo +a,iijda  = 

Eo 

^  Q(E^,a) 


0  5 

The  bracket  symbol  is  defined  by  the  following  expression 


(13) 


H 


Jx,  y  {  *  N(x)  +  ly  -;4  [N(x)  N(x  -  y)]  . 


(14) 


In  physical  terms,  Q(E  ,^2).  {E  ,^2}  represents  the  difference  (positive  or 

o  o 

negative)  between  the  number  of  absorption  processes  (upwards)  and  emission 

processes  (downwards)  involving  the  states  E^  and  E^-Q.  Clearly,  if  the 

states  E  and  E  -^2  are  in  detailed  balancing  (LTE),  the  bracket  is  zero, 
o  o 

The  intensity  of  the  local  radiation  field  can  always  be  expressed,  for  a  given 

photon  energy  J2,  in  terms  of  a  Planckian  with  a  pseudo-temperature  T.  Similarly, 

the  number -density  of  particles  in  two  states  differing  in  energy  AE  =  12  can  be 

expressed  in  terms  of  a  Maxwellian  distribution,  with  a  fixed  total  number  N  of 

o 

particles  and,  again,  a  pseudo-temperature  T'.  If  the  bracket  £^,  is  zero, 

T  =  T*.  It  will  be  shown  later  that  the  situation  where  f2^a0  for  all  12  is 

equivalent  to  LTE, 
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Returning  to  the  general  discussion,  we  write  the  steady -etate  condition  (13) 


in  the  following  form: 


for  all  E  =  const* 
o 

Mathematically,  the  steady-state  condition  can  be  visualized  in  the  following 
way:  over  the  x-y  plane,  the  function  Q(x,  y).  -^x,  y 

state  condition  (15)  then  requires  that  allong  ALL  paths  outlined  in  Fig  1,  the 


may  be  plotted*  The  steady 


integral  must  vanish. 


Let  us  finally  rewrite  Eq.  (15)  in  terms  of  x  =  E  and  z  =  E, 

'  -IX/  upper  lower 

instead  of  x  and  y,  where  x,  y  and  z  are  related  by 

X  -  z  =  y.  (16) 

The  result  reads 
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(17) 


Ob 

J  Q(x,  X-2)-  J 


The  graphic  representation  similar  to  Fig  1  is  outlined  in  Fig  2,  where  now  the 
function 


Q(x,  x-z)-  J  X,  x-z  I  =  £(x,z) 


(18) 


must  be  plotted  vs.  x  and  z. 

3.  STEADY -STATE  CONDITION  FOR  THE  RADIATION  FIELD 

Before  we  discuss  the  steady-state  condtion  for  the  particles,  Eq.  (17),  any 
further,  let  us  briefly  glance  at  the  conditions  which  a  stationary  state  imposes 
on  the  radiation  field.  So  far,  the  radiation  field  entered  our  description  only 
through  I^  as  a  parameter  in  the  bracket  function  shall  see  presently 

that  in  order  to  conserve  radiative  energy,  only  a  considerably  weaker  condtion 
must  be  satisfied,  but  that  self-consistent  solutions  to  any  given  problem  are 
subject  to  a  set  of  additional  requirements,  outlined  in  Section  6. 

There  are  two  characteristic  differences  between  the  radiation  field  and  the 
particle  distribution.  First,  the  particles  are  stationary,  whereas  the  radiation 
field  consists  of  photons  entering  and  leaving  a  certain  volume  element.  Second, 
the  number  of  particles  before  and  after  a  radiative  interaction  is  conserved  in 
steady  state,  whereas  the  number  of  photons  may  increase  or  decrease  as  a 
result  of  these  radiative  interactions.  In  fact,  a  typical  problem  in  radiation 
theory  is  the  flow  of  radiation  through  a  plasma  in  which  the  total  radiative  energy 
is  of  course  conserved,  but  in  which  the  high  energy  photons,  for  example,  are 
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broken  up  into  several  low  energy  units. 


Hence,  the  only  condition  imposed  on  the  radiation  field  under  stationary 
conditions  is  that  the  total  radiative  energy  at  each  point  is  conserved.  This 
quantity  is  found  by  multiplying  the  number  of  transitions  involving  a  state  E 

o 

by  the  appropriate  energy  and  then  integrating  over  all  states  £  from  0  to  oo, 

o 

i.  e. ,  by  forming 

OQ  oo  Oo  X 

J  dz  J  dx  •  f*(x,z)  =  jdix  jdz  •  £*(x,z)  =  0,  (19) 

0  2  0  0 

where  now 

£*(x,z)  =  (x  -  z)  Q(x,  x-z)  j  X,  X-Z^  .  (20) 


Mathematically,  the  two  representations  (19)  correspond  to  the  two  possible 
orders  of  integration  in  Cartesian  coordinates  of  the  function  f^(x,  z)  over  the 
shaded  area  in  Fig  3.  The  equality  of  integrals  of  this  type  is  shown  in  texts  on 

4 

analysis,  for  instance,  in  Jeffrey* s  book  ,  Identifying  f’J^ix,  z)  with  his  f(y),  our 
coordinates  x  and  z  with  Jeffrey* s  x  and  y,  and  letting  his  t  —¥co  T  f’*'  clearly 
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vanishes  at  infinity  on  physical  grounds  j  ,  transforms  our  expressions  (19)  into 
Jeffrey's  standard  equations • 

In  physical  terms  the  first  integral  (19)  corresponds  to  counting  first  all 
transitions  leading  into  a  fixed  lower  level  E^,  and  then  varying  (or  z  )  over 
the  whole  range,  whereas  the  second  integral  in  (19)  first  counts  all  transitions 
leading  into  a  fixed  upper  state  E^,  and  then  varies  this  upper  state  E^  (or  x  ) 
from  0  to  00. 

The  steady-state  condition  for  the  radiation  field  is  clarly  weaker  than  the 
condition  for  the  particle  distribution,  Eq.  (17),  since  it  involves  only  the 
INTEGRAL  over  expressions  of  the  type  (17).  It  will,  therefore,  in  general  not 
restrict  the  validity  of  solutions  compatible  with  (17).  However,  the  energy 
condition  (19)  is  not  redundant  in  the  sense  that  it  is  automatically  fulfilled  for  all 
solutions  to  (17),  due  to  the  additional  factor  (x-z)  under  the  integral. 

Finally,  self-consistency  of  solutions  imposes  a  set  of  conditions  developed 
in  the  context  of  the  transfer  problem  in  Section  6. 

4.  TRIVIAL  SOLUTION;  LTE 

Before  we  discuss  condition  (17)  in  more  detail,  the  existence  of  a  trivial 
solution  should  be  mentioned.  This  trivial  solution  corresponds  to  f(x,  y)  =  0 
everywhere.  On  physical  grounds,  this  means 

j  X,  x-z  I  =0  for  all  x,  z  ^  (21) 

since  a  Q  which  vanishes  everywhere  implies  no  radiative  interactions  whatsoever, 
(^1)  describes  detailed  balancing  for  all  energies  and  energy  differences. 
Intuitively  it  is  obvious  that  detailed  balancing  can  only  occur  in  LTE,  To  show 
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that  the  two  concepts  are  equivalent  under  very  general  conditions  requires,  how¬ 


ever,  a  considerable  amount  of  statistical  mechanics,  and  is  far  beyond  the  scope 
5 

of  this  paper  ,  We  shall  merely  VERIFY  that  for  radiative  transitions  between 
continuous  energy  states  detailed  balancing  and  LTE  are  equivalent. 

That  LT£  is  SUFFICIENT  for  detailed  balancing  is  obvious:  inserting  a 
Maxwellian  distribution  for  N(E)  and  a  Planck  function  for  I_  makes  the  bracket 

n 

symbol  (14)  vanish  identically,  i.e. ,  leads  to  detailed  balancing. 

Reflection  is  needed  to  show  that  LTE  is  NECESSARY  for  detailed  balancing. 
For  this  purpose,  we  write  Eq.  (14), 


}  x,y  j  =  N(x)  +  ly  ^  [K(x)  N(x-y))  x  0^ 


in  the  form 


where 


K(y)  = 


f5r 


T 


Hq  ly 


+  1 


does^not  depend  on  x.  Using  the  further  abbreviation 

G(t)  =  N(t)  /  yr^ 

we  have  the  functional  equation 

G(x-y)  =  K(y)  G(x). 

We  now  expand  K(y)  around  y  =  0, 

K(y)=  K(0)*  . 

with 

K(0)  =  1, 


Because  of  £q.  (26). 


j^dK/dy 


«  =  K' 
y  =  0  ° 


(22) 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 

(29) 
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cannot  vanish  in  general,  since  (N(x)  must  contain  the  Maxwellian  distribution  for 
which  obviously  0,  In  fact,  it  can  be  shown  that  i^  =  0,  the  normaliz¬ 

ation  integral 


oo 

J  N(E)  e”  dE  =  <  e"  >  (30) 

0 

diverges,  corresponding  to  an  infinite  temperature  (see  below). 


Next,  we  expand  G(x-y)  around  x, 

dG 

G(x-y)  =  G(x)  + - -  y  + 

dx 

For  y  0  we  find  from  Eqs.  (26),  (27)  and  (31) 

dG 

G(x) - y  =  [1  +  K'jj-y]  G(x), 

dx 

with  the  solution 

-K'  •  X 

G(x)  =  const*e  . 


or 

N(E)  dE  =  const *6  o  .  •/g’^dE  . 


(31) 


(32) 


(33) 


(34) 


£q.  (34)  represents  the  Maxwell  distribution  with 


K'jj  =  1  /  KT  . 

From  Eqs.  (23)  and  (24)  it  is  now  readily,  verified  that 


(35) 


(36) 
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(36) 


which  is  the  Planckism  distribution  with  the  same  value  (35)  for  K'  .  Hence,  we 


have  verified  that  the  condition  0  is  equivalent  to  LT£« 

5.  DISCUSSION  OF  THE  STEADY-STATE  CONDITION  FOR  PARTICLES 
We  return  now  to  the  general  steady-state  condition  for  particles,  Eq,  (17), 
with  the  additional  condition  (19)  representing  conservation  of  radiative  energy.  A 
detailed  discussion  of  sufficient  and  necessary  conditions  to  be  satisfied  by  I^ 
and  N(E)  in  order  to  be  compatible  with  Eqs.  (17)  and  (19)  will  be  given  in  a  later 
report.  It  suffices  here  to  mention  certain  MINIMUM  requirements  concerning 
the  bracket  function  (14). 

Since  the  bracket  function  must  vanish  integrated  along  ALL  paths  of  the  type 
outlined  in  Fig  2,  there  is  necessarily  at  least  one  point  on  EACH  integration  path 
for  which 


j  X,  x-z  j  =0 


(37) 


All  these  points  form  a  one -dimensional  continuum  which  from  obvious  geometrical 
considerations  falls  in  either  one  of  the  two  categories: 

1,  All  z-values,  0  <  z  <  oo,  fulfill  Eq.  (37).  The  corresponding  x-values  lie 

in  a  range  ^  with  >  0,  finite.  Similarly,  the  corresponding  y- range 

will  not,  in  general,  comprise  all  possible  values. 

2.  All  y  =  x-z  values,  0  <  y  <  oo,  fulfill  Eq.  (37).  x^, 

0<z<z;x  >0,  z  >  0. 

—  —  o  o  o 

Figs  4  and  5  illustrate  schematically  the  two  cases.  As  emphasized  before,  these 
are  minimum  requirements,  and  whether  in  an  actual  case  there  exists  one  of  the 
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outlined  minimum  solutions  will  depend  on  the  form  of  Q(E,Q). 


Fig. 4,  Case  1. 


Fig.  5.  Case  2. 


Assuming  a  solution  of  the  type  illustrated  in  Fig  5,  we  find  that  1  can  be 

O 

represented  by  a  Planckian  function  with  a  temperature  T(Q)  which  may  depend 
on  Q  =  y.  At  the  same  time,  the  ratio  N(E  +  f2)/N(E)  can  be  represented  for  all 
E -values  by  a  Maxwellian  distribution  with  the  same  T{0) -value  from  above.  Only 
if  T  is  independent  of  Q  in  the  whole  plane,  thus  {  «,z}  .  0,  v^e  come  back  to 
the  LTE-solution. 

6.  RADIATIVE  TRANSFER  AND  STEADY  STATE 
So  far,  I^  and  N(E)  were  treated  as  independent  quantities.  This  procedure, 
however,  is  in  general  not  permitted.  In  fact,  in  most  practical  cases,  I^  is  the 
sum  of  all  contributions  to  the  emission  from  the  whole  plasma  surrounding  the 
specific  point  under  discussion.  The  emission  itself,  on  the  other  hand,  dep>end8 
on  N(E). 
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There  are  numerous  cases  one  might  think  of  v^hich  range  from  the  completely 

self-consistent  problem,  in  which  is  SOLELY  due  to  the  plasma  radiation,  to 

the  opposite  extreme,  v^here  I  contains  practically  only  OUTSIDE  radiation.  The 

•• 

general  problem  is  described  by  the  transfer  equation  which  under  steady-state 
conditions  reads ^ 


and  has  the  solution 


Ola/'*’=£i  =  ®n  •  ‘a  ’ 
00  ^ 


(38) 


(39) 


The  optical  depth  is  defined  for  Eq  (39)  so  that  at  the  plasma  point  under 

study.  The  integration  is  extended  to  infinity,  i.  e.  ,  to  the  plasma  boundary.  Any 

incoming  radiation  field  at  this  boundary  can  then  be  accounted  for  by  a  properly 

matched  source  term  S  .  I  depends,  of  course,  on  direction  through 

•• 


where  s  is  an  arbitrary  spatial  coordinate, 


(40) 

is  the  absorption  coefficient. 


In  our  present  terminology,  the  emission  coefficient 


Jq(e, 


Il)>N(E)'dE 


XI  =  const 


XI 


so  that  the  dimension  of 


erg/cm  sec  sterad.  Next, 


oo 


*0.  -  N(E)]dE 


Q  =  const 


-1 


(dimension:  cm  ),  From  Eqs  (41)  and  (42),  the  source  function 

S^=  [cm"^  sec~^  sterad 

follows  immediately. 


(41) 


(42) 


(43) 
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Hence,  in  order  to  obtain  a  self-consistent  solution,  1  shuld  be  replaced  in 
Eqs  (17)  and  (19),  i.e.,  in  the  bracket  function,  by  the  appropriate  solution  (39) 
to  the  transfer  problem,  using  Eqs  (41)  -  (43)  to  represent  the  source  function. 

The  result  is  then  a  set  of  integral  equations  that  contain  solely  N(E)  as  unknown. 

It  should  be  emphasized  that  1  cannot  be  simply  replaced  by  S  Such  a 
procedure  corresponds  to  what  may  be  called  "semi-detailed^*  balancing,  implying 

®  ,  (44) 

because  of  Eq  (39)i  or 

a  Xl  ft-  (45) 


With  the  aid  of  Eqs.  (41)  -  (43),  Eq^(45)  can  be  written  as 


oo 

^  Q(E, 


i2)-N(E)-dE 


n=  const 


eo 


I  (  Q(E,a)  N—  N(E-II)  -  N(E)  ldE| 
ft.  J  lE-ft  n  = 

ft. 


const 


(46) 


or,  in  our  previous  terminology,  as 

pO 

=  const 
ft 

for  all  n  =  y.  The  path  of  integration  is  illustrated  in  Fig  6. 
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toy/€C 


This  additional  condition  will  in  general  be  compatible  with  Eq  (17)  only  for 
X,  -  0,  i.e,,  in  LTE. 
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APPENDIX:  MASER  TYPE  AMPLIFICATION 
Our  definition  (42)  allows  for  a  very  elementary  discussion  in  general  terms 
of  the  problem  of  possible  negative  absorption,  i.e, ,  maser-type  amplification 
of  continuous  radiation.  This  problem  has  recently  been  treated  by  several 
authors  in  view  of  possible  applications^. 

Since  the  factor  in  front  of  the  integral  is  positive,  the  problem  of  negative 
's  reduces  to  the  study  of 

_ , 

J  Q(E,Q)  N(E-II)  -  N(E)  ]  dE  a  R(E^)  -  R(n)  (48) 

a 

In  (481  fhe  infinite  upper  limit  is  replaced  by  a  formal  cut-off  parameter  E^  in 
order  to  ensure  finite  normalization  for  strongly  non-Maxwellian  distribution 
functions.  Negative  absorption  occurs,  if 


1r(e„)1^|r(a)|  . 


The  inequality  (49)  can,  theoretically,  be  fulfilled  in  two  cases:  since  Q  on 
physical  grounds  is  always  positive,  only  the  slope  as  a  function  of  E  is  of 
importance.  On  the  other  hand,  slope  AND  sign  must  be  taken  into  account  in 


the  case  of  the  bracket 


[]• 


In  fact,  for  a  distribution  of  the  mathematical  form 


of  the  Maxwell  distribution,  Eq  (4),  j  goes  through  zero  for  T  oo,  and  becomes 

negative  for  formally  negative  values  of  T. 

Since  all  known  interaction  mechanisms  lead  to  transition  probabilities  that 

effectively  increase  ^^dth  increasing  £,  a  combination  of 

-|n( 

(  1  <  0.  Q  K  E 

is  of  little  practical  interest.  Hence,  negative  absorption  is  confined  to  hypoth¬ 
etical  distribution  functions  whose  energy  variation  is  of  the  form 

i  +  ®  (51) 


N(E)  «(  E 


5  >  0  . 
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